LEINARTAS'S PARTIAL FRACTION DECOMPOSITION 



ALEXANDER RAICHEV 



Abstract. These notes describe Lemartas's algorithm for multivariate partial fraction 
decompositions and employ an implementation thereof in Sage. 



1. Introduction 

In [Lci78], Leinartas gave an algorithm for decomposing multivariate rational expres- 
sions into partial fractions. In these notes I re-present Lemartas's algorithm, because it 
is not well-known, because its English translation [Le!78] is difficult to find, and because 
it is useful e.g. for computing residues of multivariate rational functions; see [AY83, 
Chapter 3] and [RW12]. 

Along the way I include examples that employ an open-source implementation of 
Lemartas's algorithm that I wrote in Sage [S + 12]. The code can be downloaded from 
my website and is currently under peer review for incorporation into the Sage codebase. 

For a different type of multivariate partial fraction decomposition, one that uses iter- 
ated univariate partial fraction decompositions, see [Sto08]. 

2. Algorithm 

Henceforth let K be a field and K its algebraic closure. We will work in the factorial 
polynomial rings K[X\ and where X = Xi, . . . , Xd with d > 1. Lemartas's 

algorithm is contained in the constructive proof of the following theorem, which is [Lei78, 
Theorem 1]*. 

Theorem 2.1 (Lemartas decompositon). Let / = p/q, where p,q G iiT[X]. Let q = 

Qi 1 ' ' 'Qrn be the unique factorization of q in and let Vi = {x G K : qi(x) = 0}, 

the algebraic variety of q^ over K. 

The rational expression / can be written in the form 

f = T PA 
A 1 UeA % 

where the bi are positive integers (possibly greater than the ej), the pa are polynomials 
in K[X] (possibly zero), and the sum is taken over all subsets A C {1, . . . , m} such that 
^igaVi 7^ and {qi : i G A} is algebraically independent (and necessarily \A\ < d). 

Let us call a decomposition of the form above a Leinartas decomposition. An 
immediate consequence of the theorem is the following. 
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Corollary 2.2. Every rational expression in d variables can be represented as a sum 
of rational expressions each of whose denominators contains at most d unique irreducible 
factors. □ 

Now for a constructive proof of the theorem. It involves two steps: decomposing 
/ via the Nullstellensatz and then decomposing each resulting summand via algebraic 
dependence. We need a few lemmas. 

The following lemma is a strengthening of the weak Nullstellensatz and is proved in 
[DLLMM08, Lemma 3.2]. 

Lemma 2.3 (Nullstellensatz certificate). A finite set of polynomials {q\, . . . ,q m } C 
K[X] has no common zero in K iff there exist polynomials hi, ... , h m £ K[X] such that 



1 = y^hjqj. 



i=i 

Moreover, if K is a computable field, then there is a computable procedure to check 
whether or not the qi have a common zero in K and, if not, return the hi. □ 

Let us call a sequence of polynomials hi satisfying the equation above a Nullstellen- 
satz certificate for the q^. Note that in contrast to the usual weak Nullstellensatz, here 
the polynomials hi are in K[X] and not just in 

Some examples of computable fields are finite fields, Q, finite degree extensions of Q, 
and Q. 

Applying Lemma 2.3 we get the following lemma [LeT78, Lemma 3]. 

Lemma 2.4 (Nullstellensatz decomposition). Under the hypotheses of Theorem 2.1, 
the rational expression / can be written in the form 

/ = £ PA 



UieA* 

where the pa are polynomials in K[X] (possibly zero) and the sum is taken over all 
subsets AC {1, . . . , m} such that H ie AVi ^ 0- 

Proof. If n^Vi ^ 0, then the result holds. 

Suppose now that fl^Vj = 0. Then the polynomials q^ have no common zero in K . 
So by Lemma 2.3 

1 = Mi 1 + • • • + h m q e ™ 
for some polynomials hi in if [X]. Multiplying both sides of the equation by p/q yields 

P{h\q e i H h h m q e n 



f 



y6 m 



E phj 
ei Gi p 

i=i Qi ■■■Qi ■■■Qm 

Note that phi £ K[X\. 

Next we check each summand phi/ {ql 1 • • - q%* • • • to see whether Hj^Vj ^ 0. If so, 

then stop. If not, then apply Lemma 2.3 to qf 1 , . . . q^, . . . q^ 1 . 

Repeating this procedure until it stops yields the desired result. The procedure must 
stop, because each Vj ^ since each gj is irreducible in K[X] and hence not a unit in 
K[X\. □ 
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Let us call a decomposition of the form above a Nullstellensatz decomposition. 
Example 2.5. Consider the rational expression 

_ X 2 Y + XY 2 + XY + X + Y 
* '~ XY(XY + 1) 

in Q(X, Y). Let p denote the numerator of /. The irreducible polynomials X, Y, XY +1 G 

2 

Q[X, Y] in the denominator have no common zero in Q . So they have a Nullstellensatz 
certificate, e.g. (— Y, 0,1): 

1 = (-Y)X + (0)X + {1){XY + 1). 

Applying the algorithm in the proof of Lemma 2.4 gives us a Nullstellensatz decompo- 
sition for / in one iteration: 

pj-Y) p(l) 

7 y(iy + i) xr 
-p + v 



XY + 1 XY 

■- x - Y - 1 + x¥Tl + x + Y + l + ^f 

(after applying the division algorithm) 
1 X + Y 



Notice that 



XY + 1 XY 

f 1 1 1 



X Y XY+1 

is also a Nullstellensatz decomposition for /. So Nullstellensatz decompositions are not 
unique. 

The next lemma is a classic in computational commutative algebra; see e.g. [Kay09]. 

Lemma 2.6 (Algebraic dependence certificate). Any set S of polynomials in K[X] 
of size > d is algebraically dependent. Moreover, if K is a computable field and S is 
finite, then there is a computable procedure that checks whether or not S is algebraically 
dependent and, if so, returns an annihilating polynomial over K for 5". □ 

The next lemma is [Lci78, Lemma 1]. 

Lemma 2.7. A finite set of polynomials {q±, . . . , q m } C K[X] is algebraically dependent 
iff for all positive integers e±, . . . ,e m the set of polynomials {ql 1 , ■ ■ ■ , q%?} is algebraically 
dependent. 

Proof. A set of polynomials {qi, . . . ,q m } C K[X] is algebraically independent iff the 
m x d Jacobian matrix J(qi, . . . ,q m ) := \J^\ over the vector space K(X) d has rank 
m (by the Jacobian criterion; see e.g. [ER93]) iff for all positive integers the matrix 
(^iQi* 1 §x~) = Ji^Ti ■ ■ ■ > 9m*) over the vector space K(X) d has rank m (since we are 
just taking scalar multiples of rows) iff the set of polynomials ql 1 , . . . , g^™ is algebraically 
independent (by the Jacobian criterion). 
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Moreover, if {q 1 , . . . ,q m } is algebraically dependent, then any member of the (neces- 
sarily nonempty) elimination ideal 

C^i — Qi, ■ ■ ■ i Y m — q m , Y^ 1 — Z\, . . . , Y^™ — Z m ) k[x,y,z\ H K[Zi, . . . , Z m ], 

is an annihilating polynomial for q^ 1 , . . . , g^™. Moreover a finite basis for the elimination 
ideal can be computed using Groebner bases; see e.g. [CLO07, Chapter 3]. □ 

Applying the previous two lemmas we get our final lemma [Lei78, Lemma 2]. 

Lemma 2.8 (Algebraic dependence decomposition). Under the hypotheses of Theo- 
rem 2.1, the rational expression / can be written in the form 

J TT JH ' 

A 1 lieA Qi 

where the bi are positive integers (possibly greater than the e$), the pa are polynomials 
in K[X] (possibly zero), and the sum is taken over all subsets A C {1, . . . , m} such that 
{qi : i G A} is algebraically independent (and necessarily \A\ < d). 

Proof. If {qi, . . . , q m } is algebraically independent, then the result holds. Notice that in 
this case m < d by Lemma 2.6. 

Suppose now that {qi, . . . , q m } is algebraically dependent. Then so is {qf 1 , ■ ■ ■ , q^ 1 } 
by Lemma 2.7. Let g = Y^veS ^" ^ ^IXi; • • • >^m] be an annihilating polynomial for 
{ql 1 , ■ ■ ■ , q^ 1 }, where S C N m is the set of multi-indices such that c u ^ 0. Choose a 
multi-index a G S of smallest norm ||a|| — a± + ■ ■ ■ + a m . Then at Q := (g^ 1 , . . . , q^ 1 ) 
we have 

<?(Q) = 

v£S\{a} 
. _ ~ Z^gS\{q} 

Multiplying both sides of the last equation by p/q yields 

p _ -pc u Q u 
q~ ^ c a Q a+1 

m eiVi 
-PCu -TT Qi 



^es\{«} t=l % 



Since a has the smallest norm in S it follows that for any v £ S \ {a} there exists i such 
that CKj + 1 < z/j, so that ei(a;j + 1) < e^. So for each v G S 1 \ {a}, some polynomial 
^ei(a,+i) ^ e denominator of the right side of the last equation cancels. 

Repeating this procedure yields the desired result. □ 

Let us call a decomposition of the form above an algebraic dependence decompo- 
sition. 

Example 2.9. Consider the rational expression 

_ (X 2 Y 2 + X 2 YZ + XY 2 Z + 2XYZ + XZ 2 + Y Z 2 ) 
* '~ XYZ(XY + Z) 
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in Q(X, Y, Z). Let p denote the numerator of /. The irreducible polynomials X, Y, Z, XY+ 
Z G Q[X, Y, Z] in the denominator are four in number, which is greater than the num- 
ber of ring indeterminates, and so they are algebraically dependent. An annihilating 
polynomial for them is g(A, B, C, D) = AB + C — D. 

Applying the algorithm in the proof of Lemma 2.8 gives us an algebraic dependence 
decomposition for / in one iteration: 



c a Q Q+1 

where Q = (X, Y, Z, XY + Z) and a = (0, 0, 0, 1) 

p Q(l.l,0,0) pQ(0,0,l,0) 



£(1,1,1,2) £(1,1,1,2) 
p p 



£(0,0,1,2) £(1,1,0,2) 

P , P 



Z(XY + Z) 2 XY(XY + Z) 2 ' 

Notice that in this example the exponent 2 of the irreducible factor XY + Z in the 
denominators of the decomposition is larger than the exponent 1 of XY + Z in the 
denominator of /. Notice also that 

i i i i 



X Y Z XY+Z 

is also an algebraic dependence decomposition for /. So algebraic dependence decompo- 
sitions are not unique. 

Finally, here is Lemartas's algorithm. 

Proof of Theorem 2.1. First find the irreducible factorization of q in i^T[X]. This is a 
computable procedure if K is computable. Then decompose / via Lemma 2.4. Finally 
decompose each summand of the result via Lemma 2.8. As highlighted above, the last 
two steps are computable if K is. □ 

Example 2.10. Consider the rational expression 

2X 2 Y + AXY 2 + Y 3 — X 2 — 3XY - Y 2 



XY{X + Y){Y - 1) 



in Q(X, Y). Computing a Nullstellensatz decomposition according to the proof of Lemma 2.4 
with Nullstellensatz combination 1 = 0(X) + 1(F) + 0(X + Y) - 1(Y - 1) yields 

Y 3 + X 2 -Y 2 + X X 2 Y-2X 2 -XY 
j ~ _ + X{Y -I) + (X + Y)(Y-1) + 

-2X 3 — Y 3 — 2X 2 + Y 2 2X 2 Y -Y 3 + X 2 + 3XY + Y 2 
X(X + Y) + XY(X + Y) ' 



Computing an algebraic dependence decomposition for the last term according to the 
proof of Lemma 2.8 with annihilating polynomial g(A, B, C) = A+B—C for (X, Y, X+Y) 
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yields 



2X 2 Y -Y 3 + X 2 + 3XY + Y 2 
XY(X + Y) 

2X 2 Y -Y 3 + X 2 + 3XY + Y 2 -2X 2 Y - XY 2 - X 2 - 3XY - Y 2 
~ 1 + XY~ 2 + Y 2 (X + Y) ' 

The two equalities taken together give us a Leinartas decomposition for /. 
Notice that 

f = — 1 1 1 

J ~ X + Y + X+Y + Y - 1 

is also a Leinartas decomposition of /. So Leinartas decompositions are not unique. 

Remark 2.11. In case d — 1, Leinartas decompositions are unique once the fractions 
are written in lowest terms (and one disregards summand order). To see this, note 
that a Leinartas decomposition of a univariate rational expression f — p/q must have 
fractions all of the form Pi/q^, where q = q e \ ■ ■ • q%? is the unique factorization of q in 
This is because two or more univariate polynomials are algebraically dependent 
(by Lemma 2.6). Assume without loss of generality here that deg(p) < deg(g). It 
follows that if we have two Leinartas's decompositions of p/q, then we can write them 
in the form a\jq' + a 2 /q" = bi/q' + b 2 /q", where q = q'q" with q' and q" coprime, 
deg(ai), deg(&i) < deg(g'), and deg(a2), deg^) < deg(g"). Multiplying the equality by q 
we get a±q" + a 2 q' = b\q" + b 2 q' '■ So a\ = b\ (mod q') and a 2 = b 2 (mod q"). Thus a\ = b\ 
and a 2 = b 2 . This observation used inductively demonstrates uniqueness. 

This argument fails in case d > 2, because then a Leinartas decomposition might not 
have fractions all of the form Pi/q^- 

Remark 2.12. A rational expression already with fl^Vi 7^ and {q±, . . . ,q m } alge- 
braically independent, can not necessarily be decomposed further into partial fractions. 
For example, 

/ = — — e k(Xi, x 2 , . . . , x d ), 

with m < d can not equal a sum of rational expressions whose denominators each contain 
fewer than m of the X{. Otherwise, multiplying the equation by X\X 2 ■ ■ ■ X m would yield 



for some hi G K[X] and some nonempty subset B C {l,2,...,m}, a contradiction to 
Lemma 2.3 since {Xi : i G B} have a common zero in K , namely the zero tuple. 
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